arXiv: 1509.0405Ivl [hep-ex] 14 Sep 2015 


DPF2015-457 
September 15, 2015 


Pentaquarks and Tetraquarks at LHCb 


Sheldon StoneQ 
Department of Physics 

Syracuse University, Syracuse, N.Y. 13244, USA 


Exotic resonant structures found in and 5° decays into charmonium 
in the LHCb experiment are discussed. Examination of the J/fip system in 
yl° J/ipK~p decays shows two states each of which must be composed of 
uudcc quarks, and thus are called charmonium pentaquarks. Their masses 
are 4380 ± 8 ± 29 MeV and 4449.8 ± 1.7 ± 2.5 MeV, and their corresponding 
widths (r) are 205 ± 18 ± 86 MeV, and 39 ± 5 ± 19 MeV. The preferred 
assignments are of opposite parity, with one state having spin 3/2 and the 
other 5/2. Models of internal binding of the pentaquark states are discussed. 
Finally, another mesonic state is discussed, the Z(4430)“ that decays into qV" 
and was hrst observed by the Belle collaboration in 5° —>■ q'/L+TT” decays. 
Using a sample of approximately 25,000 signal events, LHCb determines the 
to be 1+. 


PRESENTED AT 


DPF 2015 

The Meeting of the American Physical Society 
Division of Particles and Fields 
Ann Arbor, Michigan, Angnst 4-8, 2015 


^Work supported by U.S. National Science Foundation 



1 Introduction 


In 1964 Gell-Mann [^, and separately Zweig [^, proposed that hadrons were formed from 
fundamental point like fractionally charged objects now called quarks. For most of the 
last half-century all well established baryon’s could be explained by being composed of 
three quarks and mesons a quark and an anti-quark. However, in the current decade there 
have been several observations of candidate mesonic states containing two quarks and two 
anti-quarks, called tetraquarks [^, and now, as described here, the observation of two 
pentaquark candidate baryon states [^. Such states were anticipated by Gell-Mann and 
Zweig. Predictions using theoretical mechanisms in Quantum Chromo-Dynamics (QCD) 
were made first by Jaffe in 1976 for mesons, and others for baryons in 1978 [^[^. Several 
pentaquark observations made about ten years ago were all shown to be fallacious [^. 
Thus, the recent observation of two states decaying into charmonium pentaquarks, 

found in J/ip K~p decays by the LHCb experiment is surprising. 

The decay mode was first investigated because it was suggested that it could 
contribute to the background in suppressed —?• JlipK^K~ decay that was being 
searched for and subsequently observed [^. After its discovery it was used to precisely 
measure the A° baryons lifetime [^. However, one feature of the decay that was not 
addressed was an anomalous peaking structure in the J/'ipp invariant mass spectrum, 
evident in the Dalitz plot shown in Fig. [TJ While vertical bands correspond to A* K~p 
resonances, the horizontal band can only rise from structures in the J/^p mass spectrum. 
They can also be seen in the invariant mass projections shown in Fig. 



Figure 1: Invariant mass squared oi K p versus J/V’P for candidates within ±15 MeV of the T° 
mass. 
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Figure 2: Invariant mass of (a) K p and (b) J/ijjp combinations from Jj^K p decays. 

The solid (red) curve is the expectation from phase space. The background has been subtracted. 


One may wonder if the peaking in the J/ijjp mass distribution could be caused either by 
an experimental artifact or by a conspiracy of A* amplitudes. The investigations described 
here address these questions. 

2 Analysis and results 

For this study LHCb used data corresponding to 3fb“^ of integrated iuminosity in 
7 and 8 TeV pp coiiisions. The seiection criteria are thoroughiy described in the journai 
articie [^. Here I oniy give a brief summary. Events are kept (he. triggered upon) when 
they contain a J/'ip —?• decay that is detached from the origin of the primary pp 

coiiision. 

Track combinations that form J/'ip K~p candidates are considered if the hadron 

candidates are positiveiy identified in the RICH system and have significant impact 
parameters with respect to the primary pp interaction vertex. To reduce backgrounds 
transverse momentum, pt, requirements of > 500 MeV are imposed on muons and 250 MeV 
on hadrons. Requirements on the 71° candidate include a vertex < 50 for 5 degrees 
of freedom, and a flight distance of greater than 1.5 mm. The vector from the primary 
vertex to the 71° vertex must align with the 7l° momentum so that the cosine of the angle 
between them is larger than 0.999. Candidate pApT combinations are constrained to the 
J/ip mass for subsequent use. 

Then a neural network is used to reduce backgrounds while keeping the signal 
efficiency high. The variables used are the muon identification quality, the probability 
that both hadron tracks not point at the primary pp collision vertex, the scalar sum of the 
transverse momentum (pt) of the two hadrons, and variables related to the 7l° candidate 
inclnding how well all fonr tracks form a vertex, the cosine of the angle between a vector 
from the primary vertex to the 7l° vertex and the 7l° momentum vector, flight distance. 
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and pt- 

In addition, specific backgrounds from 5° and 5° decays are vetoed. These can occur 
if the particle identification fails. We remove combinations that when interpreted as 
J/'il)K^K~ fall within ±30 MeV of the mass or when interpreted as fall 

within ±30 MeV of the mass. This requirement effectively eliminates background 
from these sources and causes only smooth changes in the detection efficiencies across the 

decay phase space. Backgrounds from decays cannot contribute significantly to 
our sample. The resulting J/ipK~p mass mass spectrum is shown in Fig. There are 
26 007±166 signal candidates containing 5.4% background within ±15 MeV (±2cr) of the 
J/'il)K~p mass peak. For subsequent analysis we constrain the J/'ipK~p four-vectors to 
give the invariant mass and the momentum vector to be aligned with the measured 


direction from the primary to the vertices 13 


In this sample specific tracking artifacts were looked for including fake tracks assembled 
from mismatched upstream and downstream segments, and multiple reconstructions of 
the same track. Having found no source of tracking artifacts we proceeded to analyze the 
decay sequences represented by the Feynman diagrams shown in Fig. This endeavor 
requires a full analysis of the amplitude for each of the two decay sequences allowing for 
their mutual interference. The amplitudes are written using six independent variables; 
one is the invariant K~p mass, rriKp-, the others are decay angles. These are shown for the 
decay sequence ^ J/ipA*; A* —>■ K~p, J/ip —>■ /r+/r“ in Fig. The A* resonances are 
modeled by Breit-Wigner amplitudes except for the yl*(1405) for which a Flatte' function 
is used 14 . All other masses, e.g. mj/^p, and decay angles can be determined from these 


six quantities. 



Figure 3; Invariant mass spectrum of J/ipK p combinations, with the total fit, signal and 
background components shown as solid (blue), solid (red) and dashed lines, respectively. 
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Figure 4: Feynman diagrams for (a) —)■ J/ipA* and (b) K decay. 


There are many A* states that can be considered, and several values of the angular 
momenta that could be present in each of their decays. Not all of these states are likely to 
be produced in our final state and not all of the allowable decay angular momenta {LS 
couplings) are likely to be present. In order to make the most general description possible 
we first use all the possible states and decay angular momenta; they are listed in Table [T] 
Then data are then fit to this model which has 146 free helicity couplings, even with 
the masses and widths of the resonant states fixed to their PDG values, done in order 
to allow the fit to converge. (Variations are considered later as part of the systematic 
uncertainties.) The results of the fit are shown in Fig. The fit gives a good description 
of the A* states as can be seen in the m-Kp spectrum but fails miserably to reproduce the 
structure in mj/^p. 

Not satisfied with using all the known A* states we tried several other different 
configurations: (i) we added all the possible S* states, (ii) we added two additional A* 
allowing their masses and widths to float in the fit and allowed spins up to 5/2 with both 
parities, and (in) we added four non-resonant components with = 1/2'^, l/2“,3/2''', 
and 3/2“. None of these fits explains the data, indeed the improvements were small. 



Figure 5: Definition of the decay angles in the A* decay chain. 
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Table 1: The yl* resonances used in the different fits. Parameters are taken from the PDG 15 
We take 5/2“ for the of the yl(2585). The number of LS couplings is also listed for both 
the “reduced” and “extended” models. To fix overall phase and magnitude conventions, which 
otherwise are arbitrary, we fix the lowest angular momentum for the yl(1520) decay. A zero entry 
means the state is excluded from the fit. 


State 

JP 

Mo (MeV) 

To (MeV) 

^ Reduced 

^ Extended 

4(1405) 

1/2- 

140541+3 

50.5 ±2.0 

3 

4 

4(1520) 

3/2- 

1519.5 ± 1.0 

15.6 ±1.0 

5 

6 

4(1600) 

1/2+ 

1600 

150 

3 

4 

4(1670) 

1/2- 

1670 

35 

3 

4 

4(1690) 

3/2- 

1690 

60 

5 

6 

4(1800) 

1/2- 

1800 

300 

4 

4 

4(1810) 

1/2+ 

1810 

150 

3 

4 

4(1820) 

5/2+ 

1820 

80 

1 

6 

4(1830) 

5/2- 

1830 

95 

1 

6 

4(1890) 

3/2+ 

1890 

100 

3 

6 

4(2100) 

7/2- 

2100 

200 

1 

6 

4(2110) 

5/2+ 

2110 

200 

1 

6 

4(2350) 

9/2+ 

2350 

150 

0 

6 

4(2585) 

? 

?«2585 

200 

0 

6 


> 2200 - 
CD 

^ 2000 F- 
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-o- 
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Figure 6: Results for (a) rriKp and (b) rriji^p for the extended A* model fit without states. 
The data are shown as (black) squares with error bars, while the (red) circles show the results of 
the fit. The error bars on the points showing the fit results are due to simulation statistics. 
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Figure 7: Definition of the decay angles in the decay sequence. 


Having failed to describe the data without a resonant state decaying into J/^|)p we 
added one. Now to do this we must write the matrix element for the decay sequence 
/1° K~ ^ J/'ipp in terms of the same decay angular variables as the previous 

decay sequence involving only A* decays. The decay angles before the appropriate rotations 
that put the decays in the same rest frames are shown in Fig.[^ The derivation of the matrix 
element in full mathematical detail is given in the arXiv article and the supplementary 
material for the Physical Review Letters publication [^. 

In each £t we minimize —21n£ where C represents the fit likelihood. The difference 
of A = — 21n£ between different amplitude models reflects the goodness of fit. For two 
models representing separate hypotheses, e.g. when discriminating between different 
values assigned to a state, the assumption of a distribution with one degree of 
freedom for A under the disfavored hypothesis allows the calculation of a lower limit 
on the significance of its rejection, i.e. the p-value (^. Therefore, it is convenient to 
express values of A as n^, where corresponds to the number of standard deviations 
in the normal distribution with the same p-value. When discriminating between models 
without and with states, ria- overestimates the p-value by a modest amount. Thus, we 
use simulations to obtain better estimates of the significance of the states. 

We perform separate fits for values of 1/2^, 3/2^ and 5/2^. The mass and width 
of the putative state are allowed to vary. The best fit prefers a 5/2''‘ state, which 
improves —21n£ by 215. Figure shows the projections for this fit. While the rriKp 
projection is well described, clear discrepancies in rrij/^p remain visible. 

The next step is to fit with two P^ states including their allowed interference. These 
fits were performed both with the reduced model and the extended model in order to 
estimate systematic uncertainties. Toy simulations are done to more accurately evaluate 
the statistical significances of the two states, resulting in 9 and 12 standard deviations, for 
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Figure 8: Results of the fit with one = 5/2+ P+ candidate, (a) Projection of the invariant 
mass of K~p combinations from J/^l;K~p candidates. The data are shown as (black) 

squares with error bars, while the (red) circles show the results of the fit; (b) the corresponding 
J/^jjp mass projection. The (blue) shaded plot shows the P+ projection, the other curves represent 
individual yl* states. 
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Figure 9: Fit projections for (a) rriKp and (b) mji^p for the reduced yl* model with two P+ states 
(see Table [^. The data are shown as solid (black) squares, while the solid (red) points show the 
results of the fit. The solid (red) histogram shows the background distribution. The (blue) open 
squares with the shaded histogram represent the Pc(4450)+ state, and the shaded histogram 
topped with (purple) filled squares represents the Pc(4380)+ state. Each yl* component is also 
shown. 
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lower mass and higher mass states, using the extended model which gives lower significances. 
The best fit projections are shown in Fig. Both m-Kp and the peaking structure in 
mj/^p are reproduced by the fit. The reduced model has 64 free parameters for the A* 
rather than 146 and allows for a much more efficient examination of the parameter space 
and, thus, is used for numerical results. The two states are found to have masses of 
4380 ± 8 ± 29 MeV and 4449.8 ± 1.7 ± 2.5 MeV, with corresponding widths of 205 ± 18 ± 86 
MeV and 39 ± 5 ± 19 MeV. (Whenever two uncertainties are quoted the first is statistical 
and the second systematic.) The fractions of the total sample due to the lower mass and 
higher mass states are (8.4 ± 0.7 ± 4.2)% and (4.1 ± 0.5 ± 1.1)%, respectively. The overall 
branching fraction has recently be determined to be 

B{Al J/tl;K-p) = (3.04 ± 0.04l°j®) x 10■^ (1) 

where the systematic uncertainty is largely due to the normalization procedure, leading to 
the product branching fractions: 

B(/lg ^ Pe(4380)+F:-p)B(Pe(4380)+ ^ J/i^p) = (2.56l}|^) x 10"^ 

P(%g ^ Pe(4450)+i^-p)P(Pe(4450)+ ^ J/i^p) = (l.25l°;^2j ^ ;Lo-^ (2) 


where all the uncertainties have been added in quadrature. 

The best fit solution has spin-parity values of (3/2“, 5/2'''). Acceptable solutions 
are also found for additional cases with opposite parity, either (3/2''', 5/2“) or (5/2''‘, 
3/2“). The five angular distributions are also well fit as can be seen in Fig. 10, 


The fit projections in different slices oi K p invariant mass are given in Fig. 11, In 


slice (a) the P/" states are not present, nor should they be as they are outside of the Dalitz 
plot boundary. In slice (d) both P/" states form a large part of the mass spectrum; there 
is also a considerable amount of negative interference between them. This can be seen 
better by examining the decay angle of the P+, 6p, the angle of the proton in J/'ipp rest 


frame with respect to the P+ direction transformed into its rest frame, shown in Fig. 12 


for the entire rtiKp range. The summed fit projections agrees very well with the angular 
distributions in the data showing that two interfering states are needed to reproduce the 
asymmetric distribution^ 


*It can be shown mathematically that the states need to be of opposite parity. 
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Figure 10: Various decay angular distributions for the fit with two P+ states. The data are 
shown as (black) squares, while the (red) circles show the results of the ht. Each ht component 
is also shown. The angles are defined in the text. 
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Figure 11: vnjj^^ in various intervals of for the fit with two states: (a) < 1.55 GeV, 

(b) 1.55 < vfiKj) < 1.70 GeV, (c) 1.70 < m^p < 2.00 GeV, and (d) m^p > 2.00 GeV. The data 
are shown as (black) squares with error bars, while the (red) circles show the results of the fit. 
The blue and purple histograms show the two states. See Fig. 


10 


for the legend. 


10 




















Figure 12: EiSciency corrected and background subtracted fit projections of the decay angular 
distributions for the two states and their sum. Values of cosdp^ near —1 are correlated with 
values of rriKp near threshold, while those near +1 are correlated with higher values. 


Systematic uncertainties are evaluated for the masses, widths and fit fractions of the 
states, and for the fit fractions of the two lightest and most significant A* states. 
Additional sources of modeling uncertainty that we have not considered may affect the 
fit fractions of the heavier A* states. The sources of systematic uncertainties are listed 
in Table They include differences between the results of the extended versus reduced 
model, varying the A* masses and widths, uncertainties in the identification requirements 
for the proton, and restricting its momentum, inclusion of a nonresonant amplitude in 
the fit, use of separate higher and lower A^ mass sidebands, alternate fits, varying the 
Blatt-Weisskopf barrier factor, d, between 1.5 and 4.5 GeV“^ in the Breit-Wigner mass 
shape-function, changing the angular momentum L by one or two units, and accounting 
for potential mis-modeling of the efficiencies. For the yl(1405) fit fraction we also added 
an uncertainty for the Flatte couplings, determined by both halving and doubling their 
ratio, and taking the maximum deviation as the uncertainty. 

The stability of the results is cross-checked by comparing the data recorded in 2011/2012, 
with the LHCb dipole magnet polarity in up/down configurations, decays, and /1° 

produced with low/high values of pp. The fitters were tested on simulated pseudoexperi¬ 
ments and no biases were found. In addition, selection requirements are varied, and the 
vetoes of and are removed and explicit models of those backgrounds added to the 
fit; all give consistent results. 

Further evidence for the resonant character of the higher mass, narrower, state is 
obtained by viewing the evolution of the complex amplitude in the Argand diagram . 
In the amplitude fits discussed above, the Pc(4450)+ is represented by a Breit-Wigner 
amplitude, where the magnitude and phase vary with mj/^p according to an approximately 
circular trajectory in the (ReA-^=, ImA'^'") plane, where A^‘^ is the mj/^p dependent 
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Table 2: Summary of systematic uncertainties on masses, widths and fit fractions, and A* fit 
fractions. A fit fraction is the ratio of the phase space integrals of the matrix element squared 
for a single resonance and for the total amplitude. The terms “low” and “high” correspond to 
the lower and higher mass states. 


Source 

Mo (MeV) 
low high 

Fo 

low 

(MeV) 

high 

low 

Fit fractions (%) 
high T(1405) T(1520) 

Extended vs. reduced 

21 

0.2 

54 

10 

3.14 

0.32 

1.37 

0.15 

masses & widths 

7 

0.7 

20 

4 

0.58 

0.37 

2.49 

2.45 

Proton ID 

2 

0.3 

1 

2 

0.27 

0.14 

0.20 

0.05 

10 < < 100 GeV 

0 

1.2 

1 

1 

0.09 

0.03 

0.31 

0.01 

Nonresonant 

3 

0.3 

34 

2 

2.35 

0.13 

3.28 

0.39 

Separate sidebands 

0 

0 

5 

0 

0.24 

0.14 

0.02 

0.03 

(3/2+, 5/2-) or (5/2+, 3/2-) 

10 

1.2 

34 

10 

0.76 

0.44 



1.5-4.5 GeV-i 

9 

0.6 

19 

3 

0.29 

0.42 

0.36 

1.91 

L^i ^ P+ (low/high)P- 

6 

0.7 

4 

8 

0.37 

0.16 



Lp+ P+ (low/high) ^ J/ipp 

4 

0.4 

31 

7 

0.63 

0.37 



hfo Al ^ 

11 

0.3 

20 

2 

0.81 

0.53 

3.34 

2.31 

b 

Efficiencies 

1 

0.4 

4 

0 

0.13 

0.02 

0.26 

0.23 

Ghange 71(1405) coupling 

0 

0 

0 

0 

0 

0 

1.90 

0 

Overall 

29 

2.5 

86 

19 

4.21 

1.05 

5.82 

3.89 


part of the Pc(4450)’'' amplitude. We perform an additional fit to the data using the 
reduced A* model, in which we represent the Pc(4450)''' amplitude as the combination 
of independent complex amplitudes at six equidistant points in the range ±r = 39 MeV 
around M = 4449.8 MeV as determined in the default fit. Real and imaginary parts of 
the amplitude are interpolated in mass between the fitted points. The resulting Argand 
diagram, shown in Fig. [T^a), is consistent with a rapid counter-clockwise change of the 
Pc( 4450)''' phase when its magnitude reaches the maximum, a behavior characteristic of 
a resonance. A similar study for the wider state is shown in Fig. [I^b); although the fit 
does show a large phase change, the amplitude values are sensitive to the details of the A* 
model and so this latter study is not conclusive. 


3 Models of pentaquark structure 


All models must explain the of the two states not just one. They also should predict 
properties of other yet to be observed states: masses, widths, J^’s. There are many 


explanations of the states. Let us start with tightly bound quarks ala Jaffe [^. Early 
work [^[7, 18 has been expanded upon recently using diquark models (T^. Here each pair 


of two quarks form a colored objects along with the lone antiquark. The three colors then 
form a colorless state, as illustrated in Fig. [Idlleft). 

Molecularly bound states, which also build on previous work po], have recently 
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Figure 13: Fitted values of the real and imaginary parts of the amplitudes for the baseline (3/2“, 
5/2+) fit for a) the Pc(4450)+ state and b) the Pc(4380)+ state, each divided into six bins 

of equal width between —F and +F shown in the Argand diagrams as connected points with 
error bars increases counterclockwise). The solid (red) curves are the predictions from 

the Breit-Wigner formula for the same mass ranges with M (F) of 4450 (39) MeV and 4380 
(205) MeV, respectively, with the phases and magnitudes at the resonance masses set to the 
average values between the two points around M. The phase convention is fixed by the yl(1520). 
Systematic uncertainties are not included. 


received much attention. Models trying to explain the states discussed here 
appeared [^, and even been disputed 22 . A molecular state configuration 
in Fig. [l4|(right). 


have already 
is illustrated 



Figure 14: (left) Illustration of a tightly bound P+ state, and (right) a molecularly bound state. 
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Other attempts at explaining the data are based on concepts of rescattering . These 
types of models were proposed to explain other resonances such as the ai(1260). These 
postdictions are made by constructing an amplitude that is consistent with the data in 
shape. They make no prediction of the magnitude of the amplitude, or its width, nor do 
they predict other final states where the phenomena could be encountered. Sometimes the 
phase motion is calculated. 

The ai(1260)''‘ saga is a good example even if it’s 51 years old, indeed as old as the 
quark model. Track measurements from a bubble chamber experiment using 3.65 GeV 
incident beam 7r+ mesons that reacted as 7r+p —>■ 7r“7r+7r+ were analyzed 24 . After 
restricting the data to have one tt+tt" mass combination consistent with the mass they 
obtained the Dalitz plot shown in Fig. [I5|[a). Removing events in the low mass pTv~ band, 
due to N* resonances, they were left with the resulting p^Tr~^ mass spectrum on the right. 
This can be explained by a higher mass state and a new af state at lower mass. 

Soon after the experimental publication, a “kinematic” (or rescattering) explanation 

I compare his amplitude shown in Fig. 16(b) with 
In the Deck diagram the beam pion 


was brought forward by Deck 25 


that of resonant production shown in Fig. 16 


a 


scatters off of a virtual pion producing a dipion pair plus an additional pion. Furthermore, 
the “Deck effect” amplitude could explain the shape of the p^7r~^ mass spectrum as shown 
in Fig. [T^ 

Over the next decade the ai was observed in different reactions and charged states 
(see for example |^). However, these were usually followed by appropriate “Deck effect” 
explanations. The situation did not become resolved until the Ui was found in t~ lepton 
decays which settled the issue around 1977. Note that there were partial wave analyses 
done that supported the scattering interpretation of the oi enhancement . The resonant 
nature of the ai proves that these analyses came to incorrect conclusions, so there was 
never a clear demonstration that Deck effect exists. 



Figure 15: (a) Dalitz plot for vr+p —> vr+p® events. The N* resonance band is indicated by 
horizontal dashed lines, (b) Histogram of the invariant mass-squared for events outside the 
N*. The peaks correspond in mass to the af and resonances. 
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b) Deck 



Figure 16: (a) Normal resonant production mechanism for the af in 7r'*"p Tr~^p^ events, 
(b) Production of a low mass enhancement via non-resonant p°7r+ scattering. (Adapted from 
Ref. 


26 


4 The Z(4430) tetraquark candidate 


The Belle collaboration in 2007 while examining decays found a relatively 

narrow peak, T ~ 45 MeV, in the mass spectrum with a mass of 4433±5 MeV 29 


This state being a charged charmonium resonance cannot be comprised of only two quarks 
and, therefore, must be a tetraquark state. 

This finding was disputed by the Babar collaboration in 2008. They wrote (SOl: “We 



Figure 17: Calculation of the mass spectrum using the “Deck effect” amplitude shown in 
Fig.j^b). 
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find that each J/^7r“ or 'tp{2S)Tr~ mass distribution is well-described by the reflection of 
the measured Kn mass and angular distribution structures. We see no significant evidence 
for a Z(4430)“ signal for any of the processes investigated.” 

Subsequently, in 2013 Belle performed a reanalysis using more data containing ?a2000 
signal events, and employing a fit to the decay amplitudes using two decay sequences 
one J/ipK*^, K*^ 7r~K~^ and 5° —?• Z~K~^, Z~ and allowing for 

interferences. [^. The result changed somewhat with the mass now being 4485±22l^® MeV 
and width 200^gl35 MeV, considerably larger than in their original paper. In addition 
they determined the to be preferentially 1+, although 0“, 1“, and 2“ could not be 
excluded. 

Using all 3 fb“^ of integrated luminosity available from LHC running in 2011 and 2012, 
the LHCb collaboration did a similar amplitude analysis with ~25,000 signal events. The 
Dalitz plot and its projections are shown in Fig. 18 


The fit projections are shown in the other plots and are in good agreement with the 




Figure 18: (top left) Dalitz plot for —)■ ip'ir decays, (top right) Distribution of 

compared with the total amplitude fits and individual component projections, (bottom left) 
Amplitude fit projections compared with data in the projection, (bottom right) Same as 

above but with the additional requirement that 1.0 < < 1-8 GeV^. 
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data only if a Z(4430)“ resonant component is included (upper right). Selecting events 
with 1.0 < < 1.8 GeV^ (lower right) shows an enhanced fraction of Z(4430)“. The 

measured mass is 4475 ± MeV and width 172 ± IdlgJ MeV are consistent with the 
Belle values. Further evidence for the resonant nature of this structure is given in the 
Argand plot, made in the same manner as for the pentaquark states, shown in Fig. 

Of course there have also been scattering models devised to explain the Z(4430) results. 
Some of these are based on the original mass determination of 4430 MeV. The average 
of the updated Belle and LHCb measurements though is 4456 MeVIt turns out that 
the sum of the masses of the T>*(2010) and Di(2420) resonances is close to 4430 MeV. 
Thus some papers considered t 
source of such rescatterings 




lat a decay such as B > T>*(2010)Di(2420)A' could be the 
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In another model B ^ D K followed by 

rescattering of the D " with the into gives rise to a peak in the mass distribution 
and a large change in phase which, however, runs clockwise in the Argand plane rather 
than counterclockwise 34 . It is interesting to note that one such calculation shows no 


rescattering effect 35 


^See the discussion on the form of the Breit-Wigner amplitude used by both experiments in Ref. 
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Figure 19: Fitted values of the Z(4430) amplitude in six bins shown in an Argand 

diagram (connected points with the error bars, with the mass increasing counterclockwise). The 
red curve is the prediction from the Breit-Wigner formula with a resonance mass (width) of 
4475 (172) MeV and magnitude scaled to intersect the bin with the largest magnitude centered 
at (4477 MeV)^. Units are arbitrary. The phase convention assumes the helicity-zero iF*^(892) 
amplitude to be real. 
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5 Conclusions 


After a half century of waiting, pentaquark states have been unmasked. Using a full 
amplitude fit to the —>■ J/ipK~p decay, the LHCb collaboration has demonstrated two 
states of opposite parties decaying into J/ipp one having a mass of 4380 ± 8 ± 29 MeV and 
a width of 205 ± 18 ± 86 MeV, while the other has a mass of 4449.8 ± 1.7 ± 2.5 MeV and 
a width of 39 ± 5 ± 19 MeV. The parities of the two states are opposite with the preferred 
spins being 3/2 for one state and 5/2 for the other. 

These states have appeared after the observation of several candidate tetraquark meson 
states. The state studied with a full amplitude analysis, the Z(4430)“, has a resonant 
amplitude with a phase change consistent with a Breit-Wigner shape as do the pentaquark 
candidates. The detailed binding mechanism of these states are subject to further studies. 
This work will lead to a better understanding of the strong interactions. Here lattice gauge 
calculations of the stability and masses of these states would be very useful. Previous 
theoretical models indicated that the presence of exotic states can modify the expected 
cooling rates of neutron stars 36 , especially for lighter mass states. Perhaps other 


implications will be revealed by further studies. 
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